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In this work we calculate the equation of state of nuclear matter for different proton fractions at 
zero and finite temperature within the Thomas Fermi approach considering three different parameter 
^SJ ' sets: the well-known NL3 and TMI and a density dependent parametrization proposed by Typel 

^— »' , and Wolter. The main differences are outlined and the consequences of imposing beta-stability in 

' these models are discussed. 
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> ■ I. INTRODUCTION AND FORMALISM 
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i Understanding the properties of nuclear matter at both, normal and high densities, is of crucial im- 
portance in explaining the appearance of neutron stars after the supernova explosion and the formation 



\ of transiron elements in nuclear reactions. 

' One of the most popular relativistic models is the non- linear Walecka model ^ , which can be used 
in order to obtain different equations of state (EOS) as far as different parameter sets are employed. In 
this work we investigate the consequences in the EOS when the parametrizations of the well known NL3 
model j^], which is a good parametrization in describing finite nuclei properties, are changed to the density 
^ • dependent one proposed in M . The new parametrization is determined by fitting several nuclear matter 
^ . bulk properties and also some finite nuclei. Both models were investigated considering two types of proton 
. fractions: fixed ones and those arising when /3-equilibrium is incorporated. Some considerations are also 
^ \ done in relation to the TMI parameter set The extension of the density dependent parametrization 
■ ■ ' to finite temperature EOS is also investigated. This extension has been partially studied in |^ for 
symmetric nuclear matter only. In nuclear collisions involving stable or radioactive neutron rich nuclei, 
in experiments yielding nuclear multifragmentation, in protoneutron stars, among inumerous examples, 
the resulting matter is known to carry a reasonable amount of isospin asymmetry. Hence, in the present 
work, a more complete and detailed study is performed in order to account for asymmetric matter as 
well. 

We start from the lagrangian density of the relativistic non-linear model, adapted in order to accom- 
modate the NL3, TMI forces and the density dependent meson- nucleon coupling constants W: 
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where cr, and pf^ are the scalar-isoscalar, vector-isoscalar and vector-isovector meson fields respectively, 
i^fj.f — df_iUJu — diyujfj,, and p^,y = 9pp,y — d^Pf^ — ^p{Pf_i x Pv)- Besides this, M is the nucleon mass, ma, TOw, 
•nip are the masses of the mesons and T^, F^^, Fp are the nucleon-meson coupling constants, k, A and ^ 
are the self-interacting scalar and vector coupling constants. In this work we investigate the differences 
arising from three parameter sets, namely NL3 Q, TMl and TW [Q. In the first two cases, Fo-, F^^ and 
Fp are the usual and gp. In the second case, the density dependent coupling constants are adjusted 

in order to reproduce some of the nuclear matter bulk properties, using the parametrization given in Q. 
Other possibilities for these parameters are also found in the literature In the TW parametrization 
the meson self-coupling constants k, A and ^ are zero. The nuclear matter bulk properties described by 
these three parameter sets are displayed in table 1. 

From the Euler-Lagrange equations we obtain the field equations of motion in the mean field appro- 
ximation for infinite matter, where the meson fields are replaced by their expectation values. In this 
approximation, the expectation value of the cr, oj and p meson fields are called and 5o respectively. 

The coupled equations read 
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where the rearrangement term Sg' is given by 



^0 



dp 



■pVo 



dTp bo dTa 



(6) 



and the scalar and baryonic densities are defined as 



Ps = {'ipip), 
p = (-07°^), 
P3 = {-ipl^Tsip). 
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In the following discussion we consider nuclear matter in the the mean-field approximation only for the 
TW parameter set. Due to translational and rotational invariance the lagrangian density reduces to 
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The conserved energy- momentum tensor can be derived in the usual fashion |7|: 



(11) 



Note that the rearrangement term included above and defined in eq.(Q) assures the energy-momentum 
conservation, i.e., d^^T^'^ = 0. From the energy-momentum tensor one easily obtains the hamiltonian 
operator: 

+V {^ml<g - \mlV^ - \mlhl\^ , (f2) 

where M* = M — Ta-4> and V is the volume of the system. In the above equation the rearrangement term 
cancels out. Notice that as a consequence, the energy density does not carry the rearrangement term 
either and can be written in the semi-classical Thomas-Fermi approximation as 
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Note that for the NL3 parametrization, the term k^q/^ + -^00/24 also appears in the energy density 
equation and for the TMI these two terms come together with — ^rf,Vo'/24. Following the notation in 
1^, the thermodynamic potential can be written as 

^ £ -TS - fXpPp - finP7i, (14) 

where S is the entropy density of a classical Fermi gas, T is the temperature, fj,p (iin) is the proton 
(neutron) chemical potential and pp and Pn are respectively the proton and neutron densities, calculated 
in such a way that p — pp + p„. We have 

(27r)3 
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where the distribution functions /i+ and /i_ for particles and anti-particles have to be derived in order 
to make the thermodynamic potential stationary for a system in equilibrium. After straightforward 
substitutions, eq.(14) becomes 

^ = 2 E / vp^+^(/^+ + /-) + ^^y^p + Y^oP3 + - ^v,' - "^hi 



+ ln(l - /,+ ) + In ( ) + ln(l - 
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For a complete demonstration of the above shown expressions obtained in a Thomas- Fermi approximation 
for the non- linear Walecka model, please refer to [p[. At this point, cq.(pq) is minimized in terms of the 
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distribution functions for fixed meson fields, i.e 



For the proton distribution function, the above calculation yields 



= z/j. (17) 
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where -E*(p) = \/p^ + M*^. Similar equations, with some sign differences are obtained for the anti- 
proton, neutron and anti-neutron distribution functions. The effective chemical potentials are then defined 
as 

p 

= - r^^vb ^feo - , 

< = ^in - T^V, + ^bo - (19) 
and the following equations for the distribution functions can be written: 

1 -Kexp[(S*(p) TM*)/r] 
In the above calculation we have used 
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and p3 = Pp — Pn- The proton fraction is defined as = PpI P- 
Within the Thomas-Fermi approach the pressure becomes 



In the NL3 model, the term — k(/)q/6 — A(/)q/24 is also present in ( ^l| ) and in the TMl model these terms 
are also accompanied by ^r^Vo'/24. It is also important to stress that the thermodynamical consistency 
which requires the equality of the pressure calculated from the thermodynamical definition and from the 
energy-momentum tensor, discussed in is also obeyed by the temperature dependent TW model. 
Another quantity of interest is the nuclear bulk symmetry energy discussed in [ p^ . It is usually defined 

as 
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with 5 — p^l p and which can be analytically rewritten as 
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where 

with kF = (I.Stt^/?)^/'^. The value and behavior of the symmetry energy at densities larger than nuclear 
saturation density are still not well established. This quantity is important in studies involving neutron 
stars and radioactive nuclei. In general, relativistic and non-relativistic models give different predictions 
for the symmetry energy. A comparison between the symmetry energies coming from the NL3 and the 
TW models is also discussed in the present work. 



II. CONSIDERING /3-STABILITY 

At this point, we introduce the ideas of fi stability and charge neutrality. In an ideal system of protons, 
neutrons, electrons and muons in equilibrium, the particle levels are filled in such a way that the /3 decays 
are forbidden. In order to study the conditions of /? equilibrium, one has to incorporate leptonic degrees 
of freedom in the lagrangian density of equation (|^), obtaining for the new lagrangian the following 
expression: 

^Ih — ^" ^leptons: C^"^) 

where 

Claptons = X! '^^ (*7m'9'^ -mi)i>i, (25) 
I 

L is given in eq. (|l|) and I describes the two lightest leptons, i.e., the electron and the muon, whose 
masses are respectively me = 0.511 MeV and = 106.55 MeV. The expressions for the energy density 
£ib and the pressure Pib are also modified by the leptons, reading: 



£ib = £ + 2Y,l -^J^+^nM + fi-) (26) 
, J (27r) 



with 



fi± = ^, J—^77p^, ^ = e,M, (27) 



l + exp[(eTw)/r] 



where fii being the chemical potentials for leptons of type I, e = ^/p^ + mf and 

Pib-P + ^T.f ^^=^(/'+ + /'-) (28) 
Stt^ ^ J ^p2 + mf 

where £ and P are given by eqs. ( |l3|) and (|2|) respectively. 

Notice that the leptons are considered as a gas of non-interacting relativistic particles, in such a way 
that the minimization of the thermodynamic potential is not altered by their presence. The already 
mentioned requirement of charge neutrality yields 



Pp = Pe+ Pti, 



(29) 
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where the electron and muon densities can be read off from equation ( |15| ) by substituting i by /. From 
the condition of chemical equilibrium in the weak processes, obtained from the minimization of the Gibbs 
potential with the conditions of baryon number and electric charge conservation, one is left with the 
following relations between the chemical potentials 

= - Me , (30) 

Mm = Me- (31) 

Common definitions for the lepton fractions are Yi — pi/p, although the lepton densities are not part of 
the baryon density. Some consequences of the imposition of f3 stability in relativistic models are discussed 




III. RESULTS AND CONCLUSIONS 

In figure |l| we show the zero temperature EOS for different proton fractions and two of the parameter 
sets used in this work, i.e., NL3 and TW. The TW parametrization makes the EOS softer not only for 
symmetric nuclear matter (Yp = 0.5), as discussed in [Q, but also for all other proton fraction possibilities. 
The same is true if (3 stability is imposed. In figure ^ the EOS is plotted for T = 10 McV and again, a 
behavior similar to that of figure |l] is observed. Notice, however, that the minima of all curves are shifted 
upwards and that the curves for Yp — 0, which do not exhibit minima for T = acquire them once the 
temperature increases. 

In figure ^ we show the EOS for neutron matter (1^ = 0) at different temperatures, namely, T = 0, 
T = 10 MeV, T = 50 MeV and T = 100 MeV. At very low densities the inchnation of the curves vary 
substancially from low to high temperatures. This is because in this region of low densities, the thermal 
energy kT is an appreciable fraction of the Fermy energy ep, making the effects of the temperature, in 
particular the particle-antiparticle creation, more dramatic in this regime than at high densities, were ep 
is greater. In figure ^ the EOS is plotted, this time for symmetric nuclear matter. One can see the change 
in the minumum from a negative to a positive value, which becomes very large for high temperatures. 
Once can also notice that the minima of all curves are slightly shifted to higer densities. 

From the analysis of figures |l|j2| we conclude that the TW parametrization is softer than the NL3 one. 
This can be explained looking for the F parametrizations in the limit of p/psat ^ 1. In this limit we have 

^i{p) — > ^^Tiipsat) ~ 0-7r^{psat) , i = <7,UJ (32) 

Ci 

Tpip) 0. (33) 

At such high densities, the system interacts mainly trough the exchange of the meson u, once the scalar 
meson a saturates as ni* — > 0. The F^^ coupling constant of the NL3 model is the same as at the 
saturation density, while equation ( ^2|) says that, in this limit, the F^^ coupling constant for the TW 
parametrization is lower than the value at the saturation point. TW is thus less repulsive at high 
densities than NL3, which makes its EOS softer. This fact has important consequences, for example 
when modeling neutron stars. A soft EOS provides a neutron star with a total mass lower than the value 
obtained with a stiff EOS. 
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We have also checked that the TW parametrization provides an EOS softer than the one obtained with 
the TMl force [|, H and closer to the relativistic Brueckner-Hartree-Fock (RBHF) EOS as can be 
seen in figures ^ and ^ for T = and T — 50 MeV and symmetric nuclear matter. The same is true 
for other proton fractions. The RBHF theory produces well the nuclear matter saturation based on the 
nucleon-nucleon interaction determined by scattering experiments. The TMl parametrization includes a 
non-linear lu term and hence works with one extra parameter which is also adjusted in order to reproduce 
nuclear matter bulk properties. We then conclude that the TW parameter set is a very useful force in 
the studies involving EOS at high densities. 

In figure ^ the symmetry energy is displayed for NL3 and TW for pure neutron matter and symmetric 
nuclear matter. Different proton fractions give rise to slightly different curves because of the difference 
in Fermi momenta and in the effective mass, which enters in E*{p). According to |To[ the symmetry 
energy at normal nuclear matter density is found to lay in between 27 and 36 MeV in the mass formula 
calculations, in the range of 28 to 38 MeV in non-relativistic models and in between 35 and 42 MeV 
in relativistic models. Notice that at the saturation point, the value for the TW parametrization (32 
MeV) is somewhat lower than for the NL3, remains in the accepted range of validity and is closer to 
the predictions of non-relativistic models. Moreover, the curves obtained for the TW model present a 
much smaller symmetry energy at larger densities and also a smoother behavior as compared with the 
curves arising from the NL3 model, which gives a more linear tendency to the curve. This result can be 
explained looking at equation (^3|), which tell us that Fp goes to zero at high densities, consequently so 
does its contribution in equation (p3[). On the opposite, F^ in the NL3 parametrization is constant as a 
function of density. 

We have finally studied the particle composition once /? stability is imposed. In figures || and ^ the 



particle composition obtained at T = respectively for NL3 and TW are shown. In figure 10 we show 
the proton and neutron composition for NL3 and TW at T = MeV and TW at T = 10 MeV. One can 
see that if the temperature does not increase much, the particle composition for a fixed parameter set 
remains basically the same. Nevertheless, it changes substancially from NL3 to TW. In particular, at 
high densities, we can see that TW is less isospin-symmetric than NL3. This is due, again, to the result 
(^), which says that, for TW, the p-nucleon interaction is suppressed at high densities. It is precisely 
this interaction which drives the systems to a more isospin-symmetric configuration at high densities, as 
we can see in figure (|^) for NL3, where this interaction survives in this limit. 

An extension of this work in order to study liquid-gas phase transition and consequent droplet formation 
is currently under investigation. 
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Table 1 Nuclear matter properties. 
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FIG. 1: Binding energy in terms of the baryon density for different proton fractions and T = Q. From top to 
bottom we show the EOS with /3-stability for NL3 (dotted hne) and for TW (long-dashed hne); Fp = for NL3 
(sohd line) and TW (large spaced dotted line); Yp = 0.3 for NL3 and TW and Yp = 0.5 for NL3 and TW. 
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FIG. 2: Binding energy in terms of the baryon density for different proton fractions and T = 10 MeV. From top 
to bottom we show the EOS with Yp = for NL3 and TW; Yp = 0.3 for NL3 and TW; Yp = 0.5 for NL3 and TW. 
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FIG. 3: Binding energy in terms of the baryon density for different temperatures, Yp = and TW. From top to 
bottom we show the EOS for T = 100 MeV, T = 50 MeV, T = 10 MeV and T = 0. 
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FIG. 4: Binding energy in terms of the baryon density for different temperatures, Yp = 0.5 and TW. From top to 
bottom we show the EOS for T = 100 MeV, T = 50 MeV, T = 10 MeV and T = 0. 
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FIG. 5: Energy density per nucleon in terms of the baxyon density for zero temperature and Yp = 0.5. From top 
to bottom we show the EOS for NL3, TMl and TW. 
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FIG. 6: Energy density per nucleon in terms of the baryon density for T = 50 MeV and Yp = 0.5. From top to 
bottom we show the EOS for NL3, TMl and TW. 
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FIG. 7: Symmetry energy in terms of the baryon density respectively for Yp = 0.0 and 0.5 for NL3 (dotted curves) 
and TW (dashed and solid curves). 
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FIG. 8: Particle composition in terms of the baryon density for T = and NL3. From top to bottom we show 
the distribution of neutrons, protons, electrons and muons. 
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FIG. 9: Particle composition in terms of the baryon density for T = and TW. Prom top to bottom we show the 
distribution of neutrons, protons, electrons and muons. 
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FIG. 10: Proton and neutron composition in terms of the baryon density. From top to bottom, looking at the 
lefthand side of the figure we show the distribution of neutrons for NL3 and T = 0, TW and T = 0, TW and 
T = 10 MeV, and for protons for TW and T = 10 MeV, TW and T = and NL3 and T = 0. 



